15 Lyapunov exponents 

Whereas fractals quantify the geometry of strange attractors, Lyaponov ex- 
ponents quantify the sensitivity to initial conditions that is, in effect, their 
most salient feature. 

In this lecture we point broadly sketch some of the mathematical issues con- 
cerning Lyaponov exponents. We also briefly describe how they are com- 
puted. We then conclude with a description of a simple model that shows 
how both fractals and Lyaponov exponents manifest themselves in a simple 
model. 

15.1 Diverging trajectories 

Lyapunov exponents measure the rate of divergence of trajectories on an 
attractor. 

Consider a flow 4>{t) in phase space, given by 




If instead of initiating the flow at 0(0), it is initiated at 0(0) +£(0), sensitivity 
to initial conditions would produce a divergent trajectory: 




Here \e\ grows with time. To the first order. 
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where 



We thus find that 



Mijit) = 



m 
d(i)j 



= M{t) e. 



(31) 



Consider the example of the Lorenz model. The Jacobian M is given by 

-P P 

M{i)= -Z{t)+r -1 -X{t) 
Y{t) X{t) -b 

We cannot solve for e because of the unknown time dependence of M{t). 
However one may numerically solve for and thus e{t), to obtain (for- 

mally) 

eit) = L{t)eiO). 



15.2 Example 1: M independent of time 

Consider a simple 3-D example in which M is time- independent. 

Assume additionally that the phase space coordinates correspond to M's 
eigenvectors. 



Then M is diagonal and 



m = 



e^i^ 
e^2t 
e^3i 



where the are the eigenvalues of M. (Recall that if f = Me, then £{t) = e^*£(0), 

where, in the coordinate system of the eigenvectors, e^* = L{t).) 

As t increases, the eigenvalue with the largest real part dominates the flow 

m- 

To express this formally, let L* be the conjugate (Hermitian) transpose of L, 
i.e. 
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Also let 

Tr(L) = diagonal sum = L^-. 

i=j 

Then 

Tr[L*(t)L(t)] = e(^^+^^)* + 6^^^+^^)* + e^^'^^*^)' 

Define 

A = Inn i In (^Tr[L*(t)L(t)]^ 

A is the largest Lyapunov exponent. Its sign is crucial: 

A < =^ £{t) decays exponentially 
A > =^ e{t) grows exponentially. 



15.3 Example 2: Time-dependent eigenvalues 



Now suppose that M{t) varies with time in such a way that only its eigen- 
values, but not its eigenvectors, vary. 



Let 



X{t) 
Y{t) 



and consider small displacements 5X{t)^5Y{t)^5Z{t) in the reference frame 
of the eigenvectors. 

Then, analogous to equation (31), and again assuming that phase space co- 
ordinates correspond to M's eigenvectors. 



5X{t) ' 




- A[m 










' 6X{t) ' 


5Y{t) 














5Y{t) 


6Z{t) _ 










cm)] . 




6Z{t) 



Here A, B, C are the time-dependent eigenvalues (assumed to be real). 
The solution for 6X{t) is 



6X{t) = 6X{0) exp 



dt'Am')] 
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1 



L 



dt'A[(P{t')] 



t 6X{0) 



t 



The RHS represents the time-average of the eigenvalue A. We assume that 
for sufficiently long times this average is equivalent to an average of A for all 
possible flows evaluated at the same time. 

In other words, we assume that the flow is ergodic. 

We denote this average by angle brackets: 



(A) is one of the three Lyapunov exponents for 4'{t). 

More sophisticated analyses show that the theory sketched above applies to 
the general case in which both eigenvectors and eigenvalues vary with time. 

15.4 Numerical evaluation 

Lyaponov exponents are almost always evaluated numerically. 

The most obvious method is the one used in the problem sets: For some e{0), 
numerically evaluate e{t) , and then flnd A such that 



0-average of 74[0(t)] 



time-average of 





£{t) \ ~ |£(0)|e 



This corresponds to the deflnition of (A) above. 
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A better method avoids saturation at the size of the attractor by successively 
averaging small changes over the same trajectory: 



8(2x) 




Here e is renormalized at each step such that 

r(r) = £(0)e^^^ 

The largest Lyaponov exponent is given by the long-time average: 

A = lim — > 7i = lim — 7 In \e{iT\ 

n. — >r)0 rt ' ^ n. — >no TIT < ^ 



n— >oo 77, 



i=l 



n^cxD nr 



Experimental data poses greater challenges, because generally we have only 
a single time series X{t). 

One way is to compare two intervals on X(^), say 

\UM\ and 
where is nearly the same on both intervals. 

Then the comparison of beyond ^2 and i!^ may yield the largest Lyaponov 
exponent. 

Another way is to reconstruct phase space by, say, the method of delays. 
Then all trajectories that pass near a certain point may be compared to see 
the rate at which they diverge. 
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15.5 Lyaponov exponents and attractors in 3-D 



Consider an attractor in a 3-D phase space. There are 3 Lyaponov exponents. 
Their signs depend on the type of attractor: 



Type 


Signs of Lyapunov exponents 


Fixed point 


(-,-,-) 


Limit cycle 


(-,-.0) 


Torus r2 


(-,0,0) 


Strange attractor 


(-0,+) 



If the attractor is a fixed point, all three exponents are negative. 

If it is a limit cycle with one frequency, only two are negative, and the third 
is zero. The zero-exponent corresponds to the direction of flow — which can 
neither be expanding nor contracting. 

Of the other cases in the table below, the most interesting is that of a strange 
attractor: 

• The largest exponent is, by definition, positive. 

• There must also be a zero-exponent corresponding to the flow direction. 

• The smallest exponent must be negative — and of greater magnitude than 
the largest, since volumes must be contracting. 



146 



15.6 Smale's horseshoe attractor 

We have seen that 

• Lyaponov exponents measure "stretching." 

• Fractal dimensions measure "folding." 

Smale's horseshoe attractor exemplifies both, and allows easy quantification. 
Start with a rectangle: 

A B 



C D 

Stretch by a factor of 2; squash by a factor of 1/(277), l > ^■ 

Now fold like a horseshoe and put back in ABCD: 
A B 




c 



D 



Now iterate the process. Stretch and squash: 
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Fold and place back in ABCD: 




Each dimension is successively scaled by its own multiplier, called a Lyaponov 
number. 

Ai = 2 {x — stretch) 
A2 = — (y - squash) 



Area contraction is given by 



The Lyapunov exponents are 



A1A2 = 1/77. 



Ai = InAi 
A2 = In A2 



Note also that vertical cuts through the attractor appear as the early itera- 
tions of a Cantor set. 

To obtain the fractal dimension, we use the definition 

£-0 ln(l/e) 

Taking the initial box height to be unity, the £, N pairs for the number N of 
segments of length e required to cover the attractor is 
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1 

1/(2^) 
l/(2ry)2 



1/(2^) 



m 



N 

T" 

2 

4 

om 



Therefore the dimension D of the Cantor set is 



In 2ri 



The dimension D' of the attractor in the plane ABCD is 

In 2 

D' = 1 + 



In 2?^' 



where we have neglected the "bend" in the horseshoe (i.e., we've assumed 
the box's width is much greater than its height. 

Note that, 

as 77 ^ 1, D' ^ 2, 
because iterates nearly fill the plane. Conversely, 

as 77 00, Z)' — ^ 1, 

meaning that the attractor is nearly squashed to a simple line. 
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